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ABSTRACT 



When a metal matrix composite is cooled down to room temperature 
from the fabrication or annealing temperature, residual stresses are induced in 
the composite due to the mismatch of the thermal expansion coefficients 
between the matrix and fiber. A method can be derived for calculating the 
internal stresses appearing in Metal Matrix composites of A1 matrix with SiC 
particles due to differences in thermal expansion coefficients. Special attention 
is paid to creep deformation in the matrix phase. The analysis shows that 
considerable internal stresses and creep deformation appear in the composites 
when subjected to cooling. 
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I. INTRODUCTION 



A. METAL MATRIX COMPOSITES(MMCs) 

Discontinuously reinforced Meted Matrix Composites(MMCs) represent a 
group of materials that combine the strength and hardness of the reinforcing phase 
with the ductility and toughness of the matrix. Aluminum alloys(Al) reinforced 
with Silicon Carbide! SiC) in particulate, platelet, or whisker form and fabricated by 
powder metallurgy methods are receiving a great deal of attention from researchers 
and engineers. 

B. THE ANALYSIS OF MMCsf Al/SICs) IN THE STABLE MEMBER 

Previous research has evaluated the dimensional stability, and the thermal 

and mechanical properties of several Al/SiC MMCs in stable— member applications 
for missile inertial guidance system. The results reveal that, although the candidate 
materials from a powder blend of SiC and A1 alloy consolidated by Vacuum Hot 
Pressing! VHP) into cylindrical billets, followed by Hot Isostatic Pressing(HIP ) for 
full densification. have better isotropic properties, these MMCs show some micro 
creeping behavior in service. The microcreep of these MMCs will affect the 
dimensional stability of stable members, and is likely to arise from two sources: (1) 
creep of the metal matrix caused by internal stresses creep conditioned by externally 
applied stresses and/or ( 2 ) phase transformations during the creep condition. A first 
step toward understanding the cause of the dimensional stability problem is to 
analyze the influence of internal stresses. The internal stresses are thermal residual 
stresses when a Metal Matrix Composite is cooled down to room temperature from 
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the fabrication or annealing temperature. These residual stresses are induced in the 
composite due to the mismatch of the thermal expansion coefficients between the 
aluminum alloy matrix and the silicon carbide particles. The model, based on 
Eshelby's for mismatch problems with simplified linear elastic material behavior for 
both particles and matrix, has been used to solve the problem of thermal residual 
stress [ref 1]. 

C. PURPOSE 

The purpose of this thesis is to calculate the thermal residual stress of Al/SiC 
using theoretical methods and then to determine what influence this thermal 
residual stress has on creep deformation. In order to do this, the thermal residual 
stresses will be determined by focusing on elastoplastic matrix and elastic inclusions 
rather than on elastic matrix and elastic inclusions. Therefore, the case where both 
phases( matrix and inclusions ) are perfectly elastic will be treated first. Then, the 
theoretical results will be compared with previously obtained experimental data for 
Al/SiC composites. 



II. LITERATURE REVIEW 



A. MMCs 

Metal matrix composites(MMCs), including eutectic composites, are becoming 
important in applications as structural components which are to be at intermediate 
and high temperatures. When MMCs are fabricated at high temperature or annealed 
at certain high temperatures, the MMCs have undesirable properties, such as low 
tensile yield and ultimate strengths. Those results are mainly due to residual 
stresses that are caused by the mismatch of the thermal expansion coefficient 
between the matrix and fiber. The residual stresses so induced have been observed 
in tungsten fiber/copper composites and in SiC whisker/6061 A1 composites and 
were based on analyzed a one— dimensional model for a continuous fiber system or 
spherical particle system. The model, based on Eshelby's equivalent inclusion 
method, has been used to solve the problem of thermal residual stressesfref 2]. 

B. ELASTIC MATRIX AND ELASTIC INCLUSIONS 

Let's consider a two— phase suspension of finite extent in which one phase is a 

matrix while the other is in the form of inclusions. It is assumed that the inclusions 

can be treated as identical ellipsoids with corresponding axes aliped. Further, the 

inclusions are assumed to be randomly distributed in the matrix in such a way that 

the suspension is homogeneous on a macroscopic scale. The task of the present 

section is to solve Eshelby's transformation problem for the case where the 

fractional volume f of the inclusions is finite. Following Eshelby [ref 2], we consider 

first an auxiliary problem, what is the resulting stress a ^ in an inclusion when every 

>j 

inclusion, having isotropic plastic moduli equal to those of the matrix, is subjected 
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to a uniform transformation strain 6 ° To find the exact solution appears to be 

» j 

impossible since no detailed information is given about the spatial distribution of 

the inclusions. All we may find in this situation are approximate average fields. The 

average constrained strain €”*' in anv inclusion may be written in the form 

» j 



e -c = e c + 6 c 

1 j 1 j 1 J 



(1) 



C • 

Here, € is the constrained strain in a typical inclusion which would be produced if 
1 j 

it alone is embedded in an infinite elastic body. This term is given, according to 
Eshelby[ref ], as 



€ C =S € T 

lj ijkl ij 



( 2 ) 



Where S is Eshelby's tensor which depends on the aspect ratios of the ellipsoidal 
ijkl 

inclusion and Poisson's ratio. (Here and in the following we use the usual summation 

1 1 

c 

convention: the range of subscripts is form 1 to 3). €’ is a contribution due to all 

> ) 

the surrounding inclusions and the presence of free boundary; this term is entirely 

elastic, and is assumed to be effectively constant in and around the typical 

1 1 

c 

inclusion. Then, it will be a good approximation to regard 6" as the average 

i j 

elastic strain in the matrix for the contribution of any one inclusion to the average 
field; outside it is negligible when there are a great number of inclusion in the 
material. For a homogeneous elastic body in static equilibrium, a volume integral of 
any component of the elastic strain associated with an internal stress should vanish 
provided that the integration is taken throughout the body. Thus 
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( i— f ) e' C + f(e~ c -e“ T ) = 0 

• j i j > j 



( 3 ) 



From (1), (2) and (3), we have 



6^ = (l-f)S 6 J + fe. T 

l j ijkl k 1 i j 



(4) 



The average internal stress in any inclusion is given by 



-I P ( ,.C c T v 

a = C ( € - 6 ) 

i] ijkl kl kl 



( 5 ) 



Where C is the elastic moduli of the matrix (and also of the inclusions in the 
ijkl 

present situation) [ref 3]. Where the inclusions have elastic moduli differing from 

those of the matrix, the same approach as used by Eshelby [ref 2] will be applicable. 

* 

The essence of this approach is terms of equivalent inclusions. Let C be the 

i j kl 

T* 

elastic moduli of the inclusions, and € . be a uniform transformation strain of every 

1 * , . . j 

inclusion. Then the transformation strain of the equivalent inclusions, € , can be 

* j 

determined by solving a set of equations 



C ( €* c -e T ) = C * ( e x -e T * ) 

ijkl v kl kl J ijkl kl kl J 



(6) 



with ( 4 ). The average stress in any inclusion is given by ( 5 ), of course. Consider 
a composite specimen which is free of stress at a references temperature To- Let the 

matrix and the inclusions have isotropic linear thermal expansion coefficients a and 

* 

a , respectively. The elastic state when the specimen is subjected to a uniform 
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temperature change, aTo=T— To, is described by the situation where the matrix 

and the inclusions undergo uniform transformation strains oaT< 5 and 

* . . , . ij 
a Respectively. Here 6 is the Kronecker delta. The resulting state of stress 
‘i ' ij 

can be analyzed by putting 

6 T+ = {a* -a) aT S (7) 

1 J lj 

The overall thermal expansion coefficients a can be expressed as a function of a, 

* * 

a , C , C , f and the aspect ratios of the inclusions, 
ijkl ijkl 

C. ELASTOPLASTIC MATRIX AND ELASTIC INCLUSIONS 

When the temperature change is great, one should expect plastic deformation, 

either in the matrix and/or in the inclusions, which relaxes the internal stress 

accumulated. An important practical case, the case where the plastic deformation 

occurs only in the matrix, will be considered in the following. Further, for 

simplicity, we assume that the plastic strain is uniform throughout the matrix and 

the plastically nondeforming inclusions can be treated as identical spheroids. Then, 

the plastic strain in the matrix, £ P, may be written in consideration of a symmetric 

U 

mode of plastic deformation and the volume constant law, as: 




( 8 ) 



Where € p is a temperature— dependent parameter. In this case, the overall strain € 

— T 

of the composite is given, with use of Equation £ =f£ , by 

ij > j 



** 

«j 



6 



( 9 ) 



e** = oaT( 5 +e p +fe T 
»j ij ij i] 

T 

Where 6 is calculated from (4) and (6) with 
i j 

G T+ = (a* — q)aT 6 — € P (10) 

1 j ‘j ij 

Now, the temperature dependence of € p will be found through an energy— balance 

consideration. Since energy dissipative process is involved in plastic deformation, 

the decrease in the elastic energy, — <® e l> corresponding to a virtual plastic 

deformation 6G p , should be balanced against the dissipation energy. If the matrix 

material has a constant flow stress a in a uniaxial deformation, the condition is 

y 

written as: 



-ffiel=(l-fV y |te p | (<T y >0) (11) 

Where E e ] = — ^fcr'^€^ 

" ij ij * 

= C|€“ + C 2 6 p (a -a) A T 

+ C 3 (a*-a) 2 AT 2 (12) 

It is noted that when the composite is elastoplatic the following condition should be 
satisfied: 



(Eel)j > (E e l)ji 



Where the subscripts I and II refer to the entirely elastic case and the elastoplastic 
case, respectively[ref 3]. 

D. CREEP 

If a stress is suddenly applied to pure metals, some solid solutions and most 
engineering alloys at a temperature near or greater than 0-5Tm(where Tm is the 
absolute melting point of the metal or alloy) deformation proceeds as shown in Fig. 
2.1. The initial application of stress causes an instantaneous elastic strain €o to 
occur. If the stress is sufficiently high an initial plastic deformation €p also occurs. 
At low temperatures, significant deformation ceases after the initial application of 
stress and an increase in stress is required to cause further deformation. At elevated 
temperatures, deformation under a constant applied load continues with time. The 
early stage of such deformation called primary creep is characterized by an initially 
high creep rate d€/dt which gradually decreases with time. Eventually a linear 
variation of creep strain accumulation with time is observed. This steady— state 
creep region is characterized by a constant minimum creep rate. Steady— state creep 
rates depend significantly on stress and temperature and are used frequently to 
compare the creep resistance among alloys. After significant deformation in 
steady— state creep, necking occurs or sufficient internal damage in the form of voids 
or cavities accumulate to reduce the cross-sectional area resulting in an increase in 
stress and creep rate. The process accelerates rapidly and failure occurs. This region 
of the creep curve is called tertiary creep. Fig. 2.1 also shows the derivative of the 
creep curve or the creep rate curve. Fig. 2.1 represents the most common creep 
strain/time behavior and Fig. 2.2 illustrates the strain rate/time behavior at 
temperatures below approximately 0.3Tm. In Fig. 2, only transient creep is 
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observed which is characterized by a continuously decreasing creep rate that 
approaches zero as the inverse of time. Such low— temperature creep behavior is 
called logarithmic creep. Fig. 2.2 also illustrates the creep rate/time behavior of 
alloys that exhibit a continuously increasing creep rate in the early stages of creep. 
Such alloys have been designated Class 1 alloys. Creep by viscous dislocation glide 
results when dislocations glide or move in slip planes under the action of an applied 
stress. These dislocations drag along solute atoms attracted to the strain fields of 
the dislocations. In order for the dislocations to move the solute atmospheres must 
diffuse in the direction of dislocation motion. At lower temperatures! ~0.3Tm) or at 
higher stress levels, creep occurs by thermally activated dislocation glide. Under 
these conditions dislocations can overcome barriers to motion without dislocation 
climb. The localized motion of short segments of small dislocations is important in 
overcoming barriers. Because of the larger contribution of stress in thermally 
activated dislocation glide a temperature and stress dependence for creep different 
from that for dislocation climb or viscous glide will be observed. fref 61 



Creep rate, de/dt Creep strain. €; 
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Tfma 




Fig. 2.2 Creep rate vs time 

f 




Fig. 2.3 Typical creep curve 
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III. THEORETICAL BACKGROUND AND DERIVATION 



A. THEORETICAL BACKGROUND 

To simulate Metal Matrix Composites! MMCs), we consider an infinite 
isotropic body containing an ellipsoidal fiber as shown in Fig. 3.1. The infinite 
body(composite) experiences a temperature T1 — TO. where T1 is assumed to be 
larger than the relaxation temperature of the matrix Tr and of course, Tr>T0. To 
make the problem simpler, the case of one fiber embedded in the infinite body is 
considered here. 




Fig. 3.1. Theoretical model 

Based on the model in Fig 2, let us consider the case of elastic 
inclusions— elastic matrix, elastoplastic matrix— elastic inclusions and creep 
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deformation: 



1. Elastic Inclusions and Matrix 



.-C 



.-C 



First let us consider the average internal constrain strain t ' . To find 

» j 



6 ' , ) we insert Equation(l) and Equationf3) into Equation(4) and find 

i L 

-0 T T 

€ =( 1 — f)S € + fe which represents the average internal constrain in anv 

i j ij kl k 1 i j 

inclusion in terms of transformation strain. Next, the average internal stress in any 

T 

inclusion can be represented by Equation(5). Here, 6 is an unknown term. 

i j 

Therefore, to find this unknown term we must consider Equation (6) and Equation 

T* 

(4). In Equation (6), * represents the inclusion, and 6 represents the uniform 

1 ^ t 

transformation strain for every inclusion. Once we have determined € , according 

— T • • 1 * 

to the Equation 6 = f£ we can determine the overall strain of the specimen, and 
ij i j . 

. . t* 

in accordance with Equation(7) we can determine the uniform transformation € 

i j 



2. Blast oplastic Matrix and Elastic Inclusions 

The case for elastoplastic matrix— elastic inclusions is different from that 

of elastic inclusion and matrix. This is due to the fact that in this case plastic 

deformation occurs only inside the matrix, so we must consider the plastic strain in 

the matrix. Therefore we can find £ ^ in accordance with Equation(S) and 

I I 

Equation( 10). From E e j = — if<r e 1 we can represent a =F(€ and 

“ ij ij ij ij 

T* n 

6 =F(G. f ). After representing them this way, 



Differentiating this with respect to € , we can obtain the value of as shown 
below: 



€ p = “2CT ( 0 “ a * aT + 2CI* 1_f K 



( 15) 



Therefore, using the value of € to find the value 6 P, the value of can be 

P ij ij 

obtained, and inserting this value into Equation(16), we are ultimately able to 
obtain the total average stress in the matrix. [ref 9] 



< ^ >n, = C; , < € >u = -f (T _I 

ij m ijkl kl Vft ij 



( 16) 



3. Creep 

Applying the value of the internal stress(above) to Equation(13). we can 
determine creep deformation. Following from this theoretical base, let us first 
calculate mathematically the thermal residual stress of Al/SiC MMCs. 

B. DERIVATION 

1. Elastic Inclusions and Elastic Matrix 

T 

First let us calculate €. . for the case of elastic inclusions and elastic matrix. 
T ^ 

In order to determine 6 , we rewrite Equation(4) in matrix form: 

i j 



C- C 1 

6 Id 


' 


622 

& 

lit 

€55 

.€66. 


=( 1-0 



Sm s,2 S,3 S,4 S,5 s 1 6 




'€,[ ' 




'€,[ ' 


S21 S22 S23 S24 S25 S26 




€22 




€22 


S31 S32 S33 S34 S35 S36 




€35 


+ f 


€ 3 i 


S41 S42 S43 S44 S45 S46 




€4$ 




S51 S52 S53 S54 S55 S56 




€5$ 




€s| 


S 6 ] Sg2 Se 3 S 64 S 65 See _ 




. 6 6$ . 




[€eU 
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Simplifying this matrix we obtain: 



r «r 

I?? 

€55 

€66 





€11 S 1 1 f + S 1 2 €i! + S 1 3 € 3 ^ 




'€,i ' 




S 2 ,€,i +e 2 U So 2 + f ) + s 2 3 € 3 ^ 




£ 2 ^ 


(1-f) 


S 3 1 € 1 1 + S 3 2 € 2 5 + €3^ S33 + f ) 
€44! S 44 + f ) 


+ f 


€4* 




€ 5 .^( S55 + f ) 

. €6§( S 66 + 0 




€d 
. eeL 



Further, simplifying this matrix we find: 



€n = €if(l— ] f)Su + €2|Si2(l— f) + €3^13(1— f) + f€|) 
€22 = e,i(l-f)S 2 , + €2^22(1 - f) + € 3 SS 23 ( 1 ~ f) + f€2$ 
€33 = €iKl-f)S 3 , + ^2$S 32 ( 1 f) + €3^33(1— f) + f €3^ 
€44 = €4$(1— f)S44 + f€4^ 

€55 = €5S(l-f)S55 + f€5^ 

€66 = €6^1-0^66 + f€e^ 



Refer to Appendix (A) for the values of Sijkl- 
Next, let us rewrite Equation(6) in matrix form. 



G 



H 




0 0 
0 0 
0 0 




0 0 
0 0 
0 0 



i/OOO 
i/OOO 
(1-*/) 0 0 0 
0 2/i 0 0 

0 0 2/i 0 

0 0 0 2/i 

1/000 
r/000 
(l-i/) 0 0 0 

0 2/i0 0 

0 0 2/i 0 

0 0 0 2/i 



e iV 




€i{ 


e 22 




6 2 * 


fa 




€ 3 | 






€ 4 ! 






€ 5 ! 


€66 j 




. € 6fi . 


r 


* 


G ld 




€l ♦ 


62 ^ 




€ 2 |* 






^33* 








€6g 




€55* 


€66 J 




J 



( 17 ) 
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Simplifying the above matrix, we have the following: 



Gj|(l— J')J+€22* / J+€33i / J=(l— J'JGGjI + t / G€ 25 +J / G€ 3 jj— ^ ( 1 — I^HGi) 

* * 

— V — ^€3^ 

J/ )J+^33 |, J =|/ G^ll + (1— t / )G€2^+* / G€3^— I'HGii — 

(l-i/)He 2 $*- t'HGaS* 

Cn^J+Ciiz^J+fasG - fc')J=i / GGii+ iXjG 2|+(1— ^)GG3S — ldlGit — 

eT^ = 0 
€ 55=0 

€66 = 0 (18) 

T 

From Equation (17) and (18), if we simplify and solve for G , we obtain: 

KGii+LG 2 ^+MG 3 N-(l-*')HGii — i'HG 2^ — J'HGsji 
NG|)+PG 2 !kQe3&=-idlG,)*— (1— t/)HG 2 £*-*'HG3!i* 

RGit+SG25+TG3^=— r'HGil — ^HG2^ 1 — j/ )HG3^ 

€ 4 $ = € 4 $ 

€s£ = € 5 ^ 

€6^ = €6^ ( 19 ) 

Now, we have 6 unknown values and 6 Equations. Therefore, applying Equation(19) 
T 

to determine G : 
ij 

c ,j- -YL - XM - Ha'(1+„) 

K 
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